Berry's phase contribution to the anomalous Hall effect of gadolinium 



^1" 
O 

o 

(N 

Of 
< 



i 

s 

I 

C 

o 
o 



> 
in 
oo 

o 
-1- 

O 

c3 



i 

c 

o 
o 



X 
S3 



S. A. BailjQ and M. B. Salamon 
Department of Physics, University of Illinois Urbana- Champaign, Urbana, IL 61801 

(Dated: February 2, 2008) 

When conduction electrons are forced to follow the local spin texture, the resulting Berry phase can 
induce an anomalous Hall effect (AHE). In gadolinium, as in double-exchange magnets, the exchange 
interaction is mediated by the conduction electrons and the AHE may therefore resemble that of 
CrC>2 and other metallic double-exchange ferromagnets. The Hall resistivity, magnetoresistance, 
and magnetization of single crystal gadolinium were measured in fields up to 30 T. Measurements 
between 2 K and 400 K are consistent with previously reported data. A scaling analysis for the Hall 
resistivity as a function of the magnetization suggests the presence of a Berry's-phase contribution 
to the anomalous Hall effect. 

PACS numbers: 72.15.Gd,75.47.Np 



I. INTRODUCTION 



While many theories account for an anomalous Hall 
effect (AHE), proportional to the magnetization of 
a material, these theories often predict effects sig- 
nificantly smaller than those found in ferromagnetic 
materialsiii2iiMi&i§i£ An even more significant deficiency 
of the conventional theories is that they predict an 
anomalous Hall resistivity that is proportional to a power 
of the longitudinal resistivity, and in the absence of a 
metal insulator transition cannot account for an AHE 
that peaks near the Curie temperature, To- Recent 
models based on a geometric, or Berry, phase have had 
great success in describing the AHE in double-exchange 
systems (e.g., manganites and chromium dioxide) and 
pyrochloresi 7 ! 8 ! 9 ! 10 ! 11 ! 12 . 13 : 14 . 15 

The anomalous Hall effect in chromium dioxide, a 
metallic double-exchange ferromagnet^ was showniS to 
agree well with the description based on geometric phase 
first suggested by Ye, et ali£ In gadolinium, as in double- 
exchange magnets, the exchange interaction among lo- 
calized (4/) core spins is mediated by the conduction 
electrons. The anomalous Hall effect may therefore re- 
semble that of Cr02 and other metallic double-exchange 
ferromagnets. Monte-Carlo simulations predict that the 
same spin-texture excitations that cause the anomalous 
Hall effect in double-exchange systems are also intrinsic 
to Heisenberg ferromagnets^ Thus it is reasonable to 
seek to explain the anomalous Hall effect in other sys- 
tems using the same theory. 

Gadolinium has a unexpectedly large anomalous Hall 
effect^ In particular, when the applied magnetic field 
is parallel to the c-axis the anomalous Hall resistivity 
peaks at p xy « —6 fifl cm just below Tc^ This makes 
it a good candidate for showing a maximum near 2/3 of 
its saturation magnetization as chromium dioxide does. 
Since gadolinium is metallic even above Tc , conventional 
theories cannot explain a maximum in the Hall effect near 
the transition temperature. 
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FIG. 1: Gadolinium resistivity vs. temperature. 



II. SAMPLE PREPARATION 

A c-axis oriented gadolinium (99.99% purity) single 
crystal was purchased from MaTecK GmbH. Two cuts 
were made parallel to an in-plane axis direction, the 
sides were polished lightly to clean up rough edges from 
the saw cuts, and the c-axis plane was thinned as much 
as possible. The resulting shape is a rectangular prism 
with an approximately square cross-section, and irregu- 
lar ends. Gold contact pads were sputtered onto the sides 
of the sample. 



III. EXPERIMENTAL RESULTS 

Data were taken using a Quantum Design Physical 
Property Measurement System (PPMS) in fields up to 
7 T. The zero field resistivity for the gadolinium crystal 
is shown in Fig.^ An alternating current (37 Hz) is ap- 
plied along the a-axis. An abrupt change in slope occurs 
at the ferromagnetic transition temperature. The resid- 
ual resistivity ratio (-R300 K./R4.2 k) is 31. For Hall effect 
and magnetization measurements, the field was applied 
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FIG. 2: (Color online) Gadolinium Hall resistivity vs. reduced 
magnetization. 



along the c-axis. The demagnetizing factor N = 0.5, and 
the saturation magnetization is 7.7 /is/Gd. The large 
values of the Hall conductivity and the magnetization al- 
lowed for very precise measurements. Fig. [21 shows the 
Hall resistivity plotted vs. reduced magnetization [m = 
M / M saturation)) these data were collected in fields up to 
7 T. The Hall resistivity increases rapidly with magneti- 
zation below the Curie temperature as domains are swept 
out. There is some indication that the data maximize at 
\Pxy\ ~ 7 fj.fl cm when |m| w 0.7. It is conventional 
to separate the Hall resistivity into ordinary (OHE) and 
anomalous (AHE) contributions: p xy = R o B in + R s fi M , 
where B m = fJ.oH apphed + /i (l - N)M. R and R s are 
the ordinary and anomalous (or spontaneous) Hall coef- 
ficients respectively. The upturns at large values of m 
are from the OHE, which is small, but not completely 
negligible. 

It is difficult to make a reliable separation of the OHE 
and AHE contributions. To obtain the values shown in 
Fig. |3| we first choose the anomalous Hall coefficient, 
R s . Next, the corresponding term (linearly proportional 
to magnetization) is subtracted from the dataset. Then, 
a linear least squares fit of Hall resistivity versus inter- 
nal field is made. The value chosen for the anomalous 
Hall coefficient is adjusted until the fitting error is mini- 
mized. The best-fit anomalous Hall coefficients are shown 
in Fig. 01 This method works even slightly above Tc , be- 
cause of the large demagnetizing correction; At temper- 
atures significantly above Tc the magnetization curves 
become linear in field, and this method fails. The other 
disadvantage to this method is that the Berry's phase 
theories predict that the anomalous Hall resistivity is 
linear in magnetization near Tc only for low values of 
rn. 

The low temperature ordinary Hall coefficient agrees 
with previously reported values (see Fig.|3J). 19 The qual- 
itative behavior is also similar. R. S. Lee and S. Legvold 
report that the ordinary Hall coefficient of gadolinium 
has temperature dependence which differs dramatically 
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FIG. 3: (Color online) Gadolinium ordinary Hall coefficient 
vs. temperature. This represents the best fit to the data in 
fields below 7 T. R. S. Lee and S. Legvold's data are shown 
for comparison^- 
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FIG. 4: (Color online) Gadolinium anomalous Hall coeffi- 
cient vs. temperature. Previously reported data are shown 
for comparison^-*— 



from those of lutetium and yttrium and cannot be ex- 
plained by a two-band model; ~ They obtained a Hall co- 
efficient which changes sign near 130 K (instead of 260 K, 
as seen in Fig. |2J and decreases even more rapidly as Tc 
is approached. The most likely cause of these discrepan- 
cies is a problem with the separation of OHE and AHE. 
Lee and Legvold only applied 3 T, whereas the values 
reported here include data up to 7 T. Indeed, when a 
subtraction was attempted using the noisier 30 T data 
(see below), the ordinary Hall coefficient did not appear 
to change sign until Tc . There are two possible explana- 
tions for this behavior. The simplest is that the AHE is 
underestimated, and the residual gives an apparent con- 
tribution to the OHE. The other possibility is that the 
sign change and the sharp increase in the magnitude of 
the Hall coefficient are real effects (possibly due to ex- 
change splitting of the conduction band). In this case, 
the decreasing magnitude that we observe at higher fields 
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FIG. 5: (Color online) Temperature scaled Hall resistivity vs. 
reduced magnetization. 



and higher magnetization may be the result of an anoma- 
lous Hall effect that is not strictly linear in magnetization 
at high fields. This non-linearity, if real, would support 
the hypothesis that Berry's phase effects contribute to 
the anomalous Hall effect in gadolinium. This contribu- 
tion would decrease as the magnetization increases, thus 
giving rise to the apparent field dependence of the ordi- 
nary Hall coefficient. The ordinary Hall resistivity will 
not be subtracted from plots of the data because of this 
dilemma. 

If the anomalous Hall effect results from the ther- 
mal excitation of topological excitations, it is possible to 
use scaling relations for the magnetization and expected 
Skyrmion density to obi a:n: |:: 
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T = p° xy T c m[l - B(x)i 



,(1-gO/01 



(1) 



0, making p xy a function of m 



where T)(x) is a scaling function of the scaling variable 
x = t/h 1 ^- 13 ^, and t and h are the reduced tempera- 
ture and magnetic field respectively. Along the critical 
isotherm t = 1 - T/T c 
only. 

In an effort to extend the results in the vicinity of 
the Curie temperature to larger values of m, we mea- 
sured both the Hall resistivity and magnetization at the 
National High Magnetic Field Laboratory in fields up 
to 30 T. The high field data are consistent with those 
taken in the PPMS, but are noisier due to problems 
both with the vibrating sample magnetometer and with 
pick-up from ripple in the Bitter magnets. Nonetheless, 
there is a clear tendency for the Hall resisitivity to reach 
an extremal value close to m = 2/3. This is shown in 
Fig. where the closed symbols in the legend are from 
the PPMS measurements and the remainder from the 
30 T experiment. The solid line is the Skyrmion expres- 
sion (Eq.^) using the critical exponents for gadolinium 21 
and D(0) = f . 



IV. DISCUSSION 

Clearly the data in Fig.[5]do not collapse well, yet sug- 
gest a tendency to fit the Skyrmion picture. The initial 
slope at Tc, p xy = — 15 pfl cm, depends on the Skyrmion 
density and spin-orbit constant through: 



rxy 



1 $q X so n e aS 
ne 7r ksTc 



(2) 



Assuming n e = 1 carrier per Gd atom, S = 7/2, and a 
Skyrmion density (n) rs 0.05 near T c , we estimate a spin- 
orbit coupling constant of X so rs 12 K. Although the fit is 
consistent with the data, the data collapse is not so good. 
The spin-orbit coupling constant also seems rather large. 
We can make a rough estimate of the spin-orbit coupling 
energy, as J. Ye et alA have done for manganite, from 
the Hamiltonian: 



so — _ o o * 



(3) 



2m 2 c 2 r dr 

Next we approximate the gradient of the potential using 



d_ v _ d Ze 2 
dr dr r 



Ze 2 

TdO. 



(4) 



where r^ is the orbital radius, and a is the lattice con- 
stant. Then an approximation of the spin-orbit coupling, 
A so , is given by: 



In the free electron model: 

k Fz = 

so 



Ze 2 h 2 k Fz 
Am 2 c 2 r<ia 



a\/3 



Ze 2 



2mc 2 rd J \2ma 2 



(5) 



(6) 



(7) 



J. Ye et al. called the middle term the "dimensionless 
coupling constant appropriate for d-orbitals," and the fi- 
nal term the "band kinetic energy."— This rough estimate 
of the spin-orbit coupling constant works out to be about 
9 K for gadolinium. 

Unlike Cr02, where only those electrons participating 
in the double-exchange contribute to the conductivity, 
Gd has both s and ci-electron contributions. It is not 
surprising, therefore, that the temperature dependence 
(below 160 K) appears to be dominated by side-jump 
processes (R s cx p 2 x )r^ as seen in a plot of R s vs. p xx in 
Fig. El A side-jump contribution, presumably from those 
portions of the Fermi surface that are not strongly spin- 
polarized, should be distinguishable from the Skyrmion 
contributions, for which R s cx e ~ E "^ kBT ^ / '(k B T)£lSL As 
a further complication, R. S. Lee and S. Legvold's data 
show a low temperature sign change of the anomalous 
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FIG. 6: (Color online) Anomalous Hall coefficient vs. resis- 
tivity squared. The residual resistivity has been subtracted. 
The ordinary Hall coefficient has been neglected when con- 
verting R. S. Lee and S. Legvold's data from R\ to R s *— 
The discrepancy in the plots is either due to an error in esti- 
mating the length between voltage contacts, or a systematic 
error in reading R. S. Lee and S. Legvold's data from their 
log-scale plot. The line is the side-jump prediction using the 
experimental coefficient for ironic- 
Hall coefficient at a temperature different from the tem- 
perature at which the anomalous Hall coefficient changes 
sign; neither side-jump nor Skyrmion models can account 
for this. Extrapolation of the contribution proportional 
to the square of the resistivity predicts a much larger 
Hall effect above 200 K than is observed. The ordinary 
Hall effect has been neglected when converting R. S. Lee 
and S. Legvold's data>i^ Berger's prediction for the side- 
jump contribution is independent of the potential, so it 
should be essentially material independent, except for the 
enhancement due to band effects^ Using the rough es- 
timate calculated for iron (see Appendix J gives a slope 
that is an order of magnitude too small for both iron and 
gadolinium^ The straight line in Fig. has a coefficient 
that is one order of magnitude larger than this estimate. 
This coefficient is consistent in magnitude with experi- 
mental values for iron between 80 K and 267 K. 1,3,5 While 
this term fits the lower temperature data, it is clearly too 
large near To 

We next explore whether the anomalous Hall effect ex- 
hibits better scaling behavior if a side-jump contribution 
is removed. We assume, arbitrarily, that a small side- 
jump process contributes to one sixth of the Hall effect 
at T c , i.e., 

RSJ = 8 

HoMq 

and that the remainder is due to Berry's-phase processes. 
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FIG. 7: (Color online) Possible Berry's phase contribution 
vs. reduced magnetization. The estimated side-jump contri- 
bution accounts for one third of the Hall effect at Tc. 



The subtracted term is strictly linear in the magnetiza- 
tion, with a temperature dependence that depends on the 
square of the zero field resistivity (excluding the residual 
resistivity). This side-jump contribution is of the same 
order of magnitude as expected theoretically, as shown 
in the appendix. Figure [7] is a plot of the difference ver- 
sus reduced magnetization, showing a better collapse of 
the data at both low and high fields, with an extremum 
in the vicinity of m = 0.6. The evidence for a decrease 
in the anomalous Hall effect at high fields is even more 
convincing after subtracting the conventional term. The 
line shown in Fig. [3 is the same as in Fig. El except the 
initial slope is reduced, and D(0) = 1 has been chosen. 
This same value for D(a;) also provided a good fit for 
Cr02^ In this case the spin-orbit coupling constant re- 
quired would be 9 K, in close agreement with our rough 
estimate. 



APPENDIX: SKEW-SCATTERING AND 
SIDE-JUMP 

More conventional explanations for the anomalous Hall 
effect include side-jump and skew scattering M> Side-jump 
scattering is when carriers scatter asymmetrically off im- 
purities. Skew scattering is a process caused by interfer- 
ence between spin-orbit coupling and second order spin- 
flip scattering. 7 In conventional ferromagnets, this theory 
yields values of R s two orders of magnitude smaller than 
experimental data (according to some authors)^ Since 
the carrier-electron spins must align with the localized 
core spins in double exchange systems, spin-flip scatter- 
ing cannot occur, and therefore skew scattering cannot 
explain the Hall effect in manganites and other systems 
with strong double exchange. 

Karplus and Luttinger developed an early model for 
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the anomalous Hall effect resulting from the spin-orbit in- 
teraction of spin-polarized conduction electrons. 23 Their 
model gave R s oc p xx , but Smit criticized their model 
arguing that a periodic potential could not cause scat- 
tering and produce the anomalous Hall effect^ Smit's 
theory, known as skew scattering, is based on anisotropic 
scattering caused by the spin-orbit interaction ; 25 After 
scattering off of an impurity, the momentum of the charge 
carriers is changed. Spin-orbit coupling makes scattering 
to one side more likely; this gives rise to the Hall effect. 
Skew scattering is generally distinguished by R s oc p xx y 
but can also give terms proportional to the square of 
the resistivity. The quadratic term occurs at high impu- 
rity concentrations (simultaneous scattering from mul- 
tiple impurities) and from phonon scattering (at least 
above the Debye temperature)^ Leribau:»2£ estimates 
the phonon scattering contribution in iron as: 

Rs = 20 ' 9 m ^T" 1 PL [l+^ 2 x 1.12 x 10-' K- 2 ] . (A.l) 

Somewhat later, Berger proposed the side-jump mech- 
anism that yields R s oc p xx £ The side-jump mechanism 
occurs when the center of mass of a carrier's wave packet 
is translated to the side while inside the scattering po- 
tential. The effect can be envisioned by picturing light 
striking a window at an angle. The refractive index of the 
window results in a displacement of the light's path, but 
no change in direction because both glass/air interfaces 
are parallel. In general, this translation can be in any 
direction, but only asymmetric (due to the spin-orbit in- 
teraction) sideways jumps will directly contribute to the 
Hall effect. Klaffky and Coleman 3 ^ estimate the side- 
jump scattering contribution in iron to be 5 times larger 



than the skew scattering contribution (Eq. IjA.ip ). an d 
given by: 

Recently, the skew and side-jump mechanisms have 
been treated simultaneously using a model based on the 
Kubo formalism and the Dirac equation^ Experimen- 
tal results for single crystal iron show that the anoma- 
lous Hall coefficient is proportional to the square of 
the resistivity between 75 K and room temperature 
The experimental coefficients, which range from 9.3 x 
10 2 f2 _1 cm _1 to 1.44 x 10 3 f2 _1 cm _1 , are much larger 
than either estimatei* 3 ^ These results do not conclu- 
sively eliminate these mechanisms as the major source of 
the anomalous Hall effect, because the estimates are only 
valid to about one order of magnitude. 
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